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' Abstract. 

pH . The Crab Nebula is hkely to be expanding into freely expanding 

\ supernova ejecta, although the energy in the ejecta may be less than is 

• typical for a Type II supernova. Pulsar nebulae much younger than the 
\ Crab have not been found and could have different properties. The search 

• for such nebulae through ultraviolet / optical line emission in core collapse 
\ supernovae, or through their X-ray emission (which could show strong 

^ ■ absorption) is warranted. Neutron stars have now been found in many 

\ young supernova remnants. There is not a clear link between neutron star 

■ and remnant type, although there is an indication that normal pulsars 
\ avoid the 0-rich remnants. In the later phases of a supernova remnant, 

■ the pulsar wind nebula is crushed by the reverse shock front. Recent 
fSI . simulations show that this process is unstable, which can lead to mixing 

' of the thermal and relativistic gases, and that the pulsar nebula is easily 

r~| . displaced from the pulsar, which can explain the position of the Vela 

Qh' pulsar relative to the Vela X radio nebula. 

o ■ 

. 

' 1. Introduction 

^ I Pulsars are thought to be born in the core collapse and explosion of massive 

. stars, so that the initial evolution and observability of the PWNe (pulsar wind 

■ " " ' nebulae) depends on the properties of the supernova and the supernova remnant 

evolution. Several phases of evolution are expected (Pacini & Salvati 1973; 
Reynolds & Chevalier 1984; Chevalier 1998). Initially the radiation from the 
PWN is absorbed by the surrounding supernova gas; high energy photons may 
be reprocessed and appear at a lower energy. After ^^lO's of years, the supernova 
becomes optically thin and the PWN radiation directly escapes. The Crab 
Nebula may be in this phase of evolution. After ~ 10^ years for a pulsar like 
that in the Crab, the spindown power from the pulsar decreases and the PWN 
fades. After ~ IQ^ years, the reverse shock wave from the supernova remnant 
moves to the center of the remnant, crushing the PWN and leading to an increase 
in brightness of the nebula. 

The discussion of PWNe has typically centered on objects like the Crab 
and its pulsar, which is estimated to have an initial spin period Pq ~ 19 ms 
and a magnetic field of 5 x 10^^ G. It is thought to belong to the main (not 
binary related) group of radio pulsars, which appear to have a fairly narrow 
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distribution of initial magnetic field strengths ~ 3 x 10^^ G (Stollman 1987). 
The initial rotation rates of these pulsars remain uncertain; studies of the pulsar 
population are influenced by selection effects (e.g., Emmering & Chevalier 1989). 
In recent years there has been increasing evidence for a class of strongly mag- 
netized neutron stars, or magnetars, with B ^ 10^^ G (Thompson & Duncan 
1995). Magnetars are thought to be associated with soft gamma-ray repeaters 
and possibly with AXPs (anomalous X-ray pulsars), both of which appear to 
contain slowly rotating neutron stars. Although the youth of AXPs is indicated 
by their association with supernova remnants, there is no evidence for PWNe. 
Finally there is a class of "quiet" neutron stars that have been found as weak 
X-ray sources in young supernova remnants, such as Gas A (Tananbaum 1999) 
and Puppis A (Petre, Becker, & Winkler 1996). Although the relative forma- 
tion rates of the normal radio pulsars, magnetars, and quiet neutron stars is not 
known, the likely youth of the objects in the latter two categories indicates that 
all three groups make a significant contribution to the total rate of massive star 
core collapses. 

Here, I discuss four topics related to neutron stars and their supernova rem- 
nants. The first is the quest for young PWNe, considerably younger than the 
Grab Nebula. Next, I consider the status of the Grab Nebula and its surround- 
ings. The third topic concerns the possible relation between a neutron star and 
its surrounding supernova, given the growing number of such associations. Fi- 
nally, I review recent investigations of the effect of a reverse shock front on a 
PWN. 



2. Searching for Young Pulsar Nebulae 

The Grab pulsar, with an age of 947 years, was until recently the youngest known 
pulsar. A newly discovered pulsar in the remnant Kes 75 has a characteristic age 
of only 723 years (Gotthelf et al. 2000) and so may be younger. However, the 
search for very young pulsars (age < 100 years) continues. The discovery of such 
objects is important for understanding the early evolution of pulsars. In their 
early phases, pulsars are expected to interact with their surrounding supernovae. 
This occurs both through the dynamical interaction of the pulsar bubble with the 
supernova and the interaction of the radiation with the surrounding supernova 
gas. If the pulsar injects energy at a constant rate E in a uniformly expanding 
central region of approximately constant density pc, the bubble radius increases 
as (Ghevaher 1977) 

where Pct^ is a constant because of the uniform expansion. The shock front 
generated by the pulsar bubble is expected to initially be radiative, with a total 
power 0.015E. The initial spin-down power of a pulsar in the vacuum dipole 
approximation can be expressed as 

Eo = 1.6 X 10» ( ^TTT^y ( ( erg s-. (2) 
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where Bp is the dipole magnetic field at the pole of the star, Rns is the neutron 
star radius, and Pq is the initial period. A dipole magnetic field perpendicular 
to the rotation axis has been assumed. The velocity of the shell at the bubble 
boundary is given by 

V = A7o(^r^ .V f , ^^f :r] ^ ( -L_X^'\m s'^ (3) 

1^1039 erg s-^y \^ 109 g s3 cm-3 y \W J ' 

where the reference value of pct^ is based on models for SN 1987A. 

The radiation can be supplemented by synchrotron radiation from the pul- 
sar bubble. The radiation from a pulsar bubble like the Crab Nebula appears 
to be from relativistic electrons and positrons that have been accelerated in a 
shock front which begins the deceleration of the pulsar wind (Rees &; Gunn 
1974; Kennel & Coroniti 1984a, 1984b). Kennel & Coroniti (1984b) find that 
the magnetization parameter a of the pulsar wind (ratio of Poynting flux to 
kinetic energy flux) must be ~ 0.003 to fit the Crab parameters. However, when 
a <^v/c flow solutions with a shock front are no longer possible; the shock wave 
moves in toward the position of the pulsar (Kennel &: Coroniti 1984a; Emmer- 
ing & Chevalier 1987). In the case of the Crab, 2000 km s-\ so that 
a ^ 6.7 X 10~3 is needed. The value indicated for the Crab is only slightly 

smaller than the limit. According to eqn. (1) with V oc t^^^, the current value 

of a was below the limit at an age of ^ 20 years, which shows that a young 
pulsar nebula may have different properties than the Crab. 

As a approaches v/c in nebular models, the fact that the shock front moves 
in to the central piilsar suggests that the rotating pulsar is in direct causal 
contact with the surrounding medium (Emmering &; Chevalier 1987; Lyutikov 
2001). There is no supersonic wind zone. The pulsar spindown law is then not 
given by eqn. (2), but is determined by the torque acting on the pulsar from the 
surrounding supernova gas through the magnetic field. The pulsar may do work 
on the surrounding medium, leading to some radiation, but the normal nebular 
synchrotron radiation does not occur. 

If a pulsar wind has a low magnetization parameter or the pulsar is powerful 
so it can drive a fast surrounding shell, it is possible that even a young pulsar 
can have a surrounding pulsar nebula with synchrotron emission from shock 
accelerated particles. Chevalier &: Fransson (1992) investigated the emission 
from a supernova that contains such a pulsar nebula. During the first 10 years, 
the emission is absorbed by the core gas aroTind the pulsar nebula and re-radiated 
at optical and ultraviolet wavelengths. The line widths depend on the pulsar 
power, but a typical width is ~ 1000 km s-^ If the core material is oxygen rich, 
the expected strong lines are [O III] AA4959,5007, O I A7774, Mg I] A4571, and 
various ultraviolet O lines. A characteristic property of the line widths is that 
they should broaden with time because as the supernova density decreases the 
ionized region becomes broader. 

The best tests for the presence of a young pulsar are in the nearest core 
collapse supernovae, SN 1987A in the LMC (Large Magellanic Cloud) and SN 
1993 J in M81. The ultraviolet/optical spectrum of SN 1987A is well explained 
by power input from 10~^ Mq of '^^Ti; any other power source is at a level 
^ 103^ erg s~^ (Chugai et al. 1997). If a quiet neutron star is present in SN 
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1987A, accretion onto the neutron star is expected, giving an Eddington limit 
luminosity. An accretion rate of at least the electron scattering Eddington rate 
is expected for ~ 200 years (Chevalier 1989), but the observed luminosity limit 
is considerably lower than the expected 2 x 10^^ erg s~^. One possibility is that a 
high Fe opacity effectively decreases the Eddington limit; a photon driven wind 
can then clear the region around a central neutron star on a timescale of 1 
yr (Fryer, Colgate, & Pinto 1999). Better limits on a central source of emission 
will become available because SN 1987A can now be spatially resolved at optical 
(Pun ct al. 2001), X-ray (Park et al. 2001a), and radio (Manchester ct al. 2001) 
wavelengths; in no case is there any sign of central emission. 

The case of SN 1993J is complicated by the fact that it is undergoing strong 
circumstellar interaction (Matheson et al. 2000). The lines that are mentioned 
above as resulting from a pulsar nebula can also be produced by circumstellar 
interaction and are observed in SN 1993J. Methods of distinguishing between the 
sources of power are that the circumstellar interaction lines tend to be initially 
broader and to narrow with time because of deceleration, whereas the pulsar 
nebula lines are initially narrower and broaden with time, as discussed above. 

While extragalactic supernovae are excellent sites to search for very young 
pulsar nebulae, somewhat older ones might be found by searching for extra- 
galactic objects like the Crab Nebula. In the LMC, the remnant 0540-69 has a 
pulsar nebula very much like the Crab and may the penultimate core collapse 
supernova in the galaxy before SN 1987A (Kirshner et al. 1989). Beyond the 
LMC, however, pulsar nebulae have not been identified. At radio wavelengths, 
Reynolds & Fix (1987) failed to find PWNe in the galaxy M33. At X-ray wave- 
lengths, deep surveys of star forming galaxies are becoming available and it is 
likely that PWNe are being detected. For example, Pence et al. (2001) have 
detected 110 sources in an image of MlOl down to a level of 10^^ erg s~^ in the 
0.5-10 keV band. The problem is to distinguish the PWNe from binary and 
other sources. A strongly absorbed spectrum is possible if absorption by the 
surrounding supernova gas is still important. 

3. The Crab and its Surroundings 

A long-standing mystery concerning the Crab Nebula is the nature of the medium 
into which it is expanding and which contains the pulsar nebula. There has not 
yet been an observation of an extended interaction region around the Crab. One 
possibility is that it is expanding into the freely expanding ejecta of SN 1054 
(e.g.. Chevalier 1977). Another is that the surroundings are a relatively slowly 
moving stellar wind from mass loss before the supernova (e.g., Nomoto et al. 
1982). 

There are a number of arguments that appear to favor the ejecta as the 
surrounding medium. In the ejecta case, the outer edge of the Crab is expected to 
be slowly accelerating or at least moving at a constant velocity. The observations 
of proper motion are consistent with free expansion, even in the extended "jet" 
on the north side of the Crab (Fesen &; Staker 1993). If the Crab were interacting 
with a slow wind, deceleration at the outer edge would be expected. 

In the ejecta interaction model, a radiative shock wave with velocity of ~ 
230 km s~^ is driven into ejecta over much of the lifetime of the Crab (Chevalier 
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&; Pransson 1992). The shocked ejecta cool and accumulate in the observed 
optical filaments. The filamentary structure can be attributed to the action of 
the Rayleigh-Taylor instability. Jun (1998) has shown in numerical simulations 
how the instability can give rise to structures like those observed in the Crab. 
Sankrit &; Hester (1997) have found observational evidence that the radiative 
shock wave is still present in the Crab, based on the strength of [O III] lines at 
the boundary of the Crab. These lines arc difficult to form by photoionization 
and are more likely to be due to a shock front. In their picture, the shock front 
is radiative only over the middle part of the Crab, so it is currently making a 
transition from the radiative to the nonradiative phase. 

A recent observation is the ultraviolet spectrum of the Crab pulsar in a 
search for absorption by freely expanding ejecta gas (Sollerman et al. 2000). An 
absorption feature in the C IV A 1550 line was found with a negative velocity 
extent to —2500 km s~^. This velocity is somewhat beyond the velocity of the 
filaments in the Crab Nebula, but not by much. It does not extend to the 
high velocities that would be needed if an explosion energy of 10^^ ergs were in 
the freely expanding ejecta with a plausible density profile. One of the initial 
arguments for expansion into free ejecta was that the ejecta could be carrying a 
normal supernova energy, as opposed to the low energy that is directly observed 
in the Crab Nebula (Chevalier 1977). 

The overall result then is that the case for interaction with freely expanding 
ejecta is good, but there is no evidence for a normal supernova energy in the 
ejecta. The lower limit found by Sollerman et al. (2000) is only ~ 1.5 x 10^^ ergs. 
Ha emission should provide a means of detecting more extended, cool ejecta, 
but only upper limits have been obtained so far (Fesen, Shull, & Hurford 1997). 
Type II supernovae with low expansion velocities have been found, such as SN 
1994W (Sollerman, Cumming, k Lundqvist 1998) and SN 1997D (Turatto et al. 
1998), but it is not clear whether these supernovae have unusually low energies. 

4. Neutron Stars and Their Supernovae 

The increasing number of neutron stars discovered in association with young su- 
pernova remnants gives the possibility of searching for correlations between core 
collapse supernovae and their central compact objects. Table 1 contains a list of 
core collapse supernovae with ages ^ 3, 000 years. The ages of 3C 58, Gil. 2-0. 3, 
and MSH 15-52 have been put in parentheses because uncertain identifications 
with SN 1181, SN 386, and SN 185, respectively, have been assumed. The first 
10 objects are in the Galaxy, the next 3 are in the LMC, and the last one is in 
the SMC. If there is not direct evidence for a compact object, the determina- 
tion of a supernova being of the core collapse variety is made based on element 
overabundances in the ejecta or in the circumstellar medium (GSM). A contro- 
versial object is Kepler's remnant, which has been suggested to have a massive 
star origin based on the circumstellar material (Bandiera 1987). Gerardy &; Fe- 
sen (2001) note the similarity between the spectrum of a circumstellar knot in 
Kepler and the quasi-stationary fiocculi in Cas A. 

The nature of the compact source and the X-ray luminosity of the PWN 
(column 3 and 4) are related in that normal pulsars appear to be capable of 
creating synchrotron wind nebulae. The quoted X-ray luminosity, Lx, is in the 
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0.2-4 keV range. The pulsar in 3C 58 was discovered only recently (Murray et 
al. 2001). The compact source in G292.0+1.8 is likely to be a normal pulsar, 
given the other correspondences with PWNe. 
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"The references to the data arc: (1) Chevahcr & Kirshner 1978; (2) Hughes et aL 2000; (3) 
Leibowitz & Danziger 1983; (4) Blanton & Helfarid 1996; (5) Gotthelf et al. 2000; (6) Fescn ct 
al. 1988 ; (7) Murray et al. 2001; (8) Davidson & Fesen 1985; (9) Murdin & Clark 1979; (10) 
Hughes ct al. 2001; (11) Torii, et al. 1997; (12) Seward et al. 1984b; (13) Tuohy & Garmire 
1980; (14) Canizares, & Winkler 1981; (15) Winkler & Kirshner 1985 ; (16) Petre et al. 1996; 
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The weak X-ray sources in Cas A, RCW 103, and Puppis A appear to be 
qualitatively different from the pulsars. In the case of Puppis A, Gaensler et 
al. (2001) have set strong limits on the presence of a radio pulsar nebula. The 
other 4 objects have neither an observed compact object nor a PWN. Three of 
these are sufficiently distant that a Cas A type source or a weak pulsar would 
not be observable. Kepler's remnant is closer and is worthy of a careful search 
for a compact object. 

The overall impression from Table 1 is that remnants with normal pulsars 
and nebulae tend not to have 0-rich ejecta, whereas the 0-rich young remnants 
tend not to have pulsars. One possible exception is 0540-69 in the LMC. This 
remnant has weak H emission, but lines of O, S, Ar, Ni, and Fc have relative 
strengths that are comparable to those observed in the Crab Nebula (Kirshner 
et al. 1989). This remnant is worthy of detailed abundance studies. A clearer 
exception is G292.0+1.8. 

The implication is that the pulsars are born in stars with a range of masses. 
The ejecta in the Crab Nebula are not 0-rich and are consistent with a low 
mass star, ~ 8 — 10 Mq (Nomoto et al. 1982). The presence of 0-rich core 
material in G292.0+1.8 implies an initial mass ^ 13 Mq. The tendency of the 
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pulsars to avoid 0-rich remnants suggests that they typically have lower mass 
progenitors, but another possibility is that the 0-rich remnants are ones with 
a dense circumstellar medium, so that the ejecta are observable through their 
interaction with the surroundings. Blair et al. (2000) suggested that N132D 
and IE 0102 had Wolf-Rayet star progenitors, and this suggestion has long 
been made for Cas A (Chevalier 1976). However, 0292.0+1.8 is a case with a 
pulsar where strong circumstellar interaction appears to be taking place (Park 
et al. 2001b). There is thus not a clear link between neutron stars and their 
supernovae, but the data are reaching the point where interesting questions can 
be asked. 

5. The Reverse Shock Wave Interaction 

In the later stage of evolution of a supernova remnant, the reverse shock front in 
the ejecta moves back to the center of the remnant, crushing a central PWN. If 
the surrounding density is ~ 1 H atom cm~^, this occurs at an age ~ 10^ years. 
Reynolds Sz Chevalier (1984) discussed the effects of the crushing on increasing 
the synchrotron luminosity of a PWN. There have recently been hydrodynamic 
studies of the crushing process. 

Van der Swaluw et al. (2001) et al. carried out one dimensional simulations 
of the process using a one fluid code. The expansion of the PWN inside the 
supernova was approximately as described by the analytical model of Reynolds 
Sz Chevalier (1984) during the pre- and post-reverse shock phases, but there were 
considerable transient waves during the transition. During the late phase, the 
radius of the PWN increases as which is close to the expansion rate of an 
adiabatic blast wave The ratio of the radius of a PWN to the radius of the 
outer shock front can thus be related to the energy injected by the pulsar, with 
only a weak dependence on age. Van der Swaluw & Wu (2001) use this property 
to obtain estimates of the total energy deposited into the PWN nebula for a 
number of observed cases and thus deduce the initial periods of the pulsars. 
The result is that the deduced periods are ^ lO's of seconds, i.e. relatively 
long. However, radiative losses from the PWN can allow a larger initial energy. 
In addition, a large energy deposited early in the evolution is primarily lost 
to adiabatic expansion, so that this method is not very sensitive to the initial 
period (Blondin, Chevalier, & Prierson 2001). 

Blondin et al. (2001) carried out two dimensional, two fluid simulations of 
the late evolution of a PWN and its surrounding supernova remnant. The PWN 
was treated as a 7 = 4/3 gas with a high sound speed. These simulations showed 
that the crushing and rc-cxpansion of the PWN are subject to hydrodynamic 
instabilities so that the thermal and nonthermal gases become mixed. The 
simulations did not include a magnetic field; the formation of magnetic filaments 
is likely to occur. In addition, the reverse shock front is unlikely to be spherically 
symmetric, so that the PWN can be pushed to the side relative to its parent 
pulsar. This effect is most likely to be observed at radio wavelengths because 
of the large ages of the radio emitting electrons, and can explain the off-center 
position of the Vela pulsar relative to the Vela X nebula. Off center radio pulsar 
nebulae have also been found in a number of other sources. 
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6. Discussion 

Of the various phases expected for pulsar nebula evolution, the earliest phase 
when nebular emission is absorbed by the surrounding supernova is least under- 
stood because of a lack of observations. The next phase when the PWN can 
be observed directly and is expanding into freely expanding supernova ejecta is 
presumably observed in many cases, although clear proof of this phase has been 
difficult to obtain even in the case of the Crab Nebula. Most of the remnants 
listed in Table 1, with ages <^ 3000 yr, are likely to be in this phase. Contin- 
ued study of these sources, especially analysis of the composition from X-ray 
emission, is warranted in order to search for correlations between neutron stars 
and their progenitor stars. In the late, post-reverse shock phase, a pulsar neb- 
ula is easily distorted and displaced by the thermal gas. Careful studies of the 
boundaries of these nebulae should reveal the nature of the interaction. 

Acknowledgments. This work was supported in part by NASA grant 
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